Correlated emission of hadrons from recombination of correlated partons 
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We discuss different sources of hadron correlations in relativistic heavy ion collisions. We show 
that correlations among partons in a quasi-thermal medium can lead to the correlated emission 
of hadrons by quark recombination and argue that this mechanism offers a plausible explanation 
for the dihadron correlations in the few GeV/c momentum range observed in Au+Au collisions at 
RHIC. 
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The recombination of thermalized quarks has recently 
been proposed as the dominant mechanism for the pro- 
duction of hadrons with transverse momenta of a few 
GeV/c in Au+Au collisions at the Relativistic Heavy 
Ion Collider (RHIC) [1 0, S 0, . While the concept 
of recombination of deconfined quarks is not new |fj , the 
RHIC data have provided compelling evidence for the 
presence of this hadronization mechanism. Quark recom- 
bination explains the enhancement of baryon emission, 
compared with meson emission, in the range of interme- 
diate transverse momenta (roughly from 2 to 5 GeV/c), 
and it provides naturally for the observed hadron species 
dependence of the elliptic flow in the same momentum 
region in terms of a universal elliptic flow curve for the 
constituent quarks Q. 

However, the model of quark recombination from a col- 
lectively flowing, deconfined thermal quark plasma ap- 
pears to be at odds with the observation of "jet-like" 
correlations of hadrons observed in the same transverse 
momentum range of 2 to 5 GeV/c ja, @- Trigger- 
ing on a hadron, e.g., with transverse momentum 2.5 
GeV/c < pt < 4 GeV/c, the data shows an enhancement 
of hadron emission in a narrow angular cone around the 
direction of the trigger hadron in a momentum window 
below 2.5 GeV/c. Can such correlations be reconciled 
with the claim that hadrons in this momentum range are 
mostly created by recombination of quarks? 

Obviously, the observation is incompatible with any 
model which assumes that no correlations exist among 
the quarks before recombination. Such correlations re- 
quire deviations from a global thermal equilibrium in the 
quark phase. One mechanism is already well established: 
a strong, but anisotropic and locally varying, collective 
flow produces correlations among hadrons after recombi- 
nation. Indeed, the hadronic elliptic flow correlation is 
known to be larger than the elliptic flow of the quarks 
before recombination 0, 0], because the parameter 
characterizing the magnitude of elliptic flow for a hadron 
is proportional to the number of its constituent quarks. 

One would generally expect that correlations among 
the quarks, when present before their recombination into 
hadrons, will be amplified by the hadronization process. 
Here, we argue that a certain amount of two-body cor- 



relations among the partons of a quark-gluon plasma 
produced in a relativistic heavy-ion collision is to be 
expected. Energetic partons produced in such a colli- 
sion lose a significant amount of energy through colli- 
sions with thermal partons on their way out of the dense 
medium, but do not completely thermalize before they 
form hadrons. The dissipated energy and momentum 
are absorbed by the surrounding medium, increasing its 
temperature slightly and setting it into motion in the 
direction of the energetic parton. This "wake effect" 
produces correlations among medium partons (S) and 
the originally energetic parton (H). If the parton loses 
enough energy to become indistinguishable from the ther- 
mal medium, all that remains is a narrowly directed, 
"jetty" flow pattern within the thermal medium. Such 
soft partons can either recombine with each other upon 
hadronization (denoted by SS for a meson) , or a soft par- 
ton can recombine with a hard jet parton (SH) [il llfj|. 
or a jet parton can fragment outside the medium (F) to 
form a hadron. 

Recently, the STAR collaboration presented experi- 
mental evidence that jet cones are not isolated from 
the surrounding medium, but that jets and the medium 
strongly influence each other They found that 

hadrons correlated with the jet participate in the lon- 
gitudinal expansion of the medium. We will not attempt 
here to give a quantitative description of how the corre- 
lations in the parton phase arise. Rather we will show 
how such correlations between partons can translate into 
hadron correlations in a hadronization scenario with re- 
combination and fragmentation. We will focus on dime- 
son production by three processes: the recombination of 
four thermal partons into two mesons (SS-SS), jet cor- 
relations by double fragmentation (F-F) and fragmenta- 
tion accompanied by a soft-hard recombination between 
a jet fragment and a thermal parton (F-SH). These pro- 
cesses are schematically depicted in Fig. Obviously, 
there are three additional possibilities for producing two 
mesons in this recombination-fragmentation picture, SH- 
SH, SH-SS and F-SS, and even more if baryons are in- 
volved. However they do not involve any essentially new 
aspects and we will not consider them further here. 

We start by discussing recombination from a thermal, 
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FIG. 1: Schematic pictures of the processes F-F, F-SH 
and SS-SS for producing two mesons A, B. Thick lines are 
medium quarks. Possible correlations are indicated by curly 
lines. 



but correlated ensemble of quarks (SS-SS), following the 
formalism described byr Fries et al. 2\. The creation 
of a pair of mesons A, B from four partons can be ex- 
pressed as a convolution of the meson Wigner functions 
&a,&b and the Wigner function W1234 of the four par- 
tons, integrated over the hadronization hypersurface E: 



E A E B d 6 N AB 
d?P A cPP B 



= C A b / d<j A da B ^A®^B®W 12 zi, (1) 



where Cab = CaCb is a degeneracy factor. The gener- 
alization to the case where A, B, or both are baryons 
is straightforward. In the absence of correlations the 
n-parton Wigner function W\... n is approximated by a 
product of classical one-parton phase space distributions 
Wi . This approximation yields satisfactory results for in- 
clusive hadron spectra. However, it should not come as 
a surprise that it is not sufficient to describe hadron cor- 
relations. For this reason, we go one step further and 
include two-parton correlations in the form 



Here tfio.n and yo,7r are the widths of the Gaussians in 
azimuth and rapidity, respectively. The two terms of the 
sum correspond to correlations initiated by an energetic 
parton (</> = 0) and its recoil partner ((p — it). Cq and c„ 
give the strength of the near side and far side correlations, 
while the functions fo,Tr(j>Ti,PTj) describe the transverse 
momentum dependence of the correlations. The func- 
tions Sb,7r(0i, <Jj) parametrize the spatial localization of 
the parton correlations on the hypersurface E. As indi- 
cated above, we assume that So,n = 1 , if o"i , <Jj G V c and 
So.n = otherwise. 

In this letter we restrict the discussion to near side cor- 
relations (cjr = 0). We also refrain from exploring the Pt 
dependence, because of the present lack of data, which 
would allow us to constrain the function fo(pn,PTj)- 
Since we work with fixed Pt windows here, for which 
experimental data exist, we shall simply set /o = 1, ab- 
sorbing all numerical factors into the parameter Co- We 
also assume that co < 1, allowing us to neglect terms of 
higher power in the correlations, such as Cq, C0V2 or v\. 

Following |3] we integrate over the spatial coordinates 
(assuming Vhadron < 14 < Vj) and the quark momenta 
transverse to the hadron momentum. Note that, un- 
like for single inclusive hadron spectra, the width of the 
hadron wave function in the azimuthal direction could, 
in principle, interfere with the correlation width <pQ. For 
simplicity, we neglect such effects here. We also use the 
narrow wave function approximation (see 0), which was 
shown to provide a good description of the measured 
spectra and elliptic flow. For the thermal parton dis- 
tributions Wi we use Boltzmann distributions with tem- 
perature T, radial flow rapidity t\t and a boost invariant 
hadronization hypersurface E at fixed proper time r Q- 

The dimeson spectrum for the SS-SS process is 



Wl234 « WlW2W3W4{l + Cij) 
i<j 



(2) 



Here Cij is the correlation function between partons i 
and j. We shall not make special assumptions about 
the origin of the correlations among partons, except that 
they do not rapidly vary in momentum, can be localized 
around a certain direction and are confined to a subvol- 
ume V c of the fireball. This is compatible with corre- 
lations being induced by the interaction of an energetic 
parton with the thermalized medium, as discussed above 
as a likely source for partonic correlations. 

As an example we are going to evaluate the azimuthal 
correlations between mesons emitted around midrapidity 
in Au+Au collisions at RHIC. We assume that momen- 
tum correlations of partons are restricted to a cone in 
rapidity y and azimuthal angle <p. Hence we choose the 
correlation functions to be of the form 
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where we introduced the abbreviation 
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with 



I + rrij are the quark masses and 

Mt is the transverse mass of the meson. This quantity is 
related to the single inclusive meson spectrum at midra- 
pidity which, including elliptic flow and two-parton cor- 
relations, is given by [2| hi(PT)(l + co)(l + 2w 2 cos(20)). 
The effect of the space-time correlation volume V c has 
been approximated in by a rescaling of the normal- 
ization constant cq = cqV c /(tAt). 
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We note that the amplification factor Q = 4 in front 
of the Gaussian term counts the number of possible cor- 
relations of a quark in meson A with a quark in meson 
B. Likewise the term 2<5o accounts for the correlations of 
quarks inside the same meson. It is then easy to see that 
the amplification factor is Q = 6 for a meson-baryon pair 
and Q = 9 for a baryon-baryon pair. This result con- 
firms our expectation that correlations within the par- 
tonic medium are amplified in the recombination process, 
similar to the amplification of the elliptic flow v 2 by the 
number of valence quarks. 

The experiments at RHIC measure the associated par- 
ticle yield per trigger hadron A. After subtracting the 
uncorrelated background and using the notation A<f> = 
\4>A — 0b|j the relevant observable is defined as 

-l ( dN A B d(N A N B y 



Y AB (A<t>) = N A > 



(6) 



\d(A(f>) d{A<j)) 

The particle yields are integrated over the kinematic win- 
dows of the hadrons A and B with exception of the rela- 
tive angle A<f>. The Pt spectra of associated particles 
have recently been studied in |l2| in a recombination 
picture. Neglecting quadratic terms of correlation co- 
efficients, the background subtraction cancels the term 
proportional to 1 + 2cq in Eq. J3J and leads to the result 

N A Y AB (A<f) = Qdoe-Wt*® N A N B /(2n), (7) 

were Ni — 2ir J dyidPTiPTihi(PTi) is the total particle 
number for species i in the kinematic window. The effect 
of a possible correlation in rapidity — not discussed here 
— can be absorbed into the normalization cq. The uncor- 
related background up to first order in the coefficients c<j, 
v 2 is given by (27r) _1 N A N B [l + 2c + 2v 2A v 2B cos(2A^)], 
where v 2 i is the average elliptic flow of hadron species i 
in the kinematic window. 

Dihadron production through fragmentation from a jet 
(F-F) is described by dihadron fragmentation functions. 
These have recently been discussed by Majumder and 
Wang [laj . We assume that they can be factorized into 
single hadron fragmentation functions D a /h with an ap- 
propriate scaling of the momentum variable. Since the 
formalism is strictly collinear, we introduce a Gaussian 
smearing in relative azimuthal angle and rapidity of the 
two hadron momenta. Integrating over rapidities we ob- 
tain for hadrons emitted around midrapidity 
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Here g a (p) — E p d 3 N a /d 3 p is the invariant spectrum of 
the fragmenting parton a. The factor / contains the in- 
tegration over the correlated hadron rapidities y A and 
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FIG. 2: YX°B e which is Yab integrated over < < 0.94, 
for meson (left panel) and baryon triggers (right panel) as a 
function of centrality. Insert: Associated yield as a function 
of A<j> before background subtraction using SS-SS and F-F 
contributions at an impact parameter 6 = 8 fhi. 



y B in their respective windows around midrapidity, as- 
suming that the parton spectra g a are slowly varying. 
Of course, the width for this mechanism does not 
generally coincide with that introduced for correlations 
among the thermal partons. The kinematic limits are 
z a = 2P TA /^/s and z\ = Pta/(Pta + Ptb)- A-E denotes 
the average energy loss of parton a before fragmentation. 

For the F-SH process we adopt the following simple 
model for dimeson production. Suppose a hard parton 
a with initial momentum p a dresses itself with a pair bb 
(where b is any flavor) and that the pair ab hadronizes 
into a meson A. Instead of further fragmentation the re- 
maining parton b can pick up a soft parton c from the 
medium and recombine into another meson B. It is clear 
that this is only the simplest possible realization of F-SH. 
We assume that the production of A can be described by 
a fragmentation function D a / A and we approximate the 
correlated emission by a product ansatz, where one of 
the parton distributions entering recombination is com- 
ing from a jet. Eventually we arrive at 

d(N A Y AB (A$)) T ,8C B M TB 

= 2-nlv 9a(Pa) 



dP TA dP TB 
x (27T^)-i/2 e -(A0) 2 /(20^) /o 
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Ptb sinh r/x 



2T 
Pta 



Pta + Ptb 12 



(9) 



with p a = Pta + Ptb /2 + AE. The normalization con- 
stant V < 1 arises from restricting the integration over £ 
to a subspace given by the jet cone. Our notation im- 
plies that the fragmented hadron A is the trigger hadron, 
but the process where B is the trigger has to be taken 
into account as well. Details of these calculations will be 
provided in a forthcoming publication. 

Fig. |21 shows numerical examples of the different con- 
tributions. We use windows of 1.7 GeV/c < Ptb < 2.5 



4 



GeV/c for associated particles and 2.5 GeV/c < Pta< 
4.0 GeV/c for trigger particles and \y\ < 0.35 as in 
Our calculation includes charged pions and kaons as well 
as protons and antiprotons. The parameters of the ther- 
mal parton phase are taken from [2j. We use the minijct 
distributions from [l4| and KKP fragmentation functions 
|l5|. The adjustable parameters in our model are the az- 
imuthal width cf>o, which is chosen to be 0.2 for all pro- 
cesses, the correlation strength Co, and v. 

Fig. compares the background subtracted associated 
yield Yab, integrated over < A(f> < 0.94 as a func- 
tion of centrality for the following scenarios: (i) F-F pro- 
cess only, (ii) FF-FF and SS-SS with no soft correlations 
(co = 0), and (hi) the same with cq — 0.08. The left 
panel is for meson triggers, the right panel for baryon 
triggers. The insert shows the uncorrelated background 
due to the associated yield from the SS-SS process and 
independent fragmentation as well as the signal for a near 
side correlation from SS-SS and F-F for c = 0.08. 

As seen in the figure, the F-F mechanism (squares) 
produces strong near-side correlations, which are larger 
for baryon triggers. However, the trigger yields from 
this process are small in the considered window, thus 
adding trigger particles from recombination dilutes the 
signal dramatically (full diamonds). This effect is even 
more pronounced for baryons (open diamonds). Switch- 
ing on SS-SS correlations strongly increases the hadron 
correlations (triangles). Note that a constant value of 
Co corresponds to a correlation volume V c scaling with 
-'Vpart, which is not likely realistic. We also show the 
result for a more realistic fixed V c corresponding to 
c = 0.08 x 100/iVp ar t (circles), which makes Yab vary 
weakly with centrality. 

We also show F-SH correlations for 7r — 7r and v = 0.5 
(stars). The yields in this case are so small that the 
contribution is negligible compared with SS-SS and F- 
F. It was already pointed out in 2] that with the parton 
parametrization found by fitting the single inclusive spec- 
tra, soft-hard processes are subdominant. Other groups 
using different parametrizations have come to different 
conclusions 0,O|. 

The PHENIX preliminary results do not show clearly 
identifiable near-side jet cones above background for 
baryon triggers J9j , which suggests that one must be cau- 
tious with the interpretation of the A0-integrated associ- 
ated yields in this case. In our model the relative behav- 
ior of baryon and meson triggers depends on the relative 
strength of the SS-SS and F-F contributions. Double 
fragmentation (F-F) would predict a larger associated 
yield for baryons triggers than for meson triggers. For 
the most realistic case (F-F and SS-SS with c « 0.08 
and fixed correlation volume) meson and baryon triggers 
result in dihadron correlations of approximately equal 
magnitude. More experimental information is needed, 
including measurements of p — p collisions, to determine 
which scenario describes the data best. 



In summary, we have shown that correlations among 
partons in a quark-gluon plasma naturally translate into 
correlations between hadrons formed by recombination 
of quarks with an amplification factor similar to the one 
obtained for elliptic flow. Preliminary data from the 
PHENIX collaboration are consistent with hadron pro- 
duction by quark recombination if two-parton correla- 
tions of order 10% within a fixed correlation volume are 
assumed to be present in the deconfined phase. We con- 
clude that the existence of localized angular correlations 
among hadrons are not in contradiction to the recom- 
bination scenario, since a plausible mechanism for the 
formation of such parton correlations — passage of hard 
partons through the dense matter — exists in nuclear 
collisions at RHIC. 
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